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Abstract 

We present explicit formulae for Weber-Schafheitlin type integrals and give them an interpretation as 
the kernel of a physically relevant operator related to the hamiltonian of Aharanov and Bohm. In particular, 
we derive explicit formulae for Weber-Schafheitlin type integrals with exponent larger or equal 1, which 
are distributions on R+. We discuss several special cases. 



1 Introduction 

Our aim is to calculate the integral 



k'^ J^(xk)Jjj{k)Ak, (1-1) 







for suitable ^,1^ £ C and x G M_|_ (0, 00), where is the Bessel function of the first kind of order /i. In 
the literature, it is known as the Weber-Schafheitlin discontinuous integral with exponent p (or —p, depending 
on the convention). 

In the case when Re p < 1 (and under some additional assumptions on p,, v) it is convergent to a function, in 
general not continuous in a; = 1. It has been derived in several ways and analysed in many special cases, for 
which we refer in particular to [ W] and [DF] . It has been applied in numerous problems, let us mention here 
only two recent works — [HN], [SS]. 

The case Re p > 1 is more delicate and requires a distributional approach. Nevertheless, it is quite natural to 
consider it; indeed, it appears in some problems where it plays an important role ([KR], [KeR]). In addition 
to that, we show in this paper that it is the kernel of an operator physically relevant for the Aharanov-Bohm 
system. There have already been successful attempts to derive useful expressions for the distributional case 
of (1.1) for special values of parameters, by Kellendonk and Richard [KR] (for p = 1), by Miroshin [M] 
(p ^ 1 asymptotic) and by Salamon and Walter [SW] (recurrence formulae and special values of parameters 
with p E Z). Motivated by the wish of exhausting all unsolved cases, we provide explicit formulae for (1.1) 
for arbitrary p with positive real part. We discuss also some special cases and compare them with the results 
mentioned earher. 

The paper is constructed as follows. Section 2 serves as an additional motivation, linking the hamiltonian 
of Aharanov and Bohm to the integral we discuss further In Section 3, we derive formulae for integrals 
involving the modified Bessel function of the second type i^T^, as in the approach of Dixon and Ferrar [DF]. 
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We then use them to compute the integrals involving the Hankel functions of the first and second kind — 

oo 

l<i''H^{xK)J^{K)dK, (1.2) 

closely related to (1.1). Our treatment of the integrals (1.2) is a generalization of the approach adopted in 
[KR]. We consider separately the cases Rep < and Rep > 0. The second one is examinated in Section 
4 and includes in particular the distributional cases. The main result for the integral (1.1) is contained in 
Proposition 1.1, and follows as a simple consequence of the computations of (1.2). 



1.1 The main result 

In the following proposition, we give formulae for the Weber-SchaftheitUn integral (1.1) with exponent 

Rc /? > 0. The result is a distribution on IR+ (in general not regular), and its form depends on whether 
p is an integer number (which is not the case if Rep < 0). We use the notation for the rescaled Gauss 
hypergeometric function 

2FI (a, 6; c; z) := y(^c+1) ^' ^' ' 

Proposition 1.1 For any p.,iy G C and Re p > satisfying Re (p + i/ + 1) > and x G M+, the integral 

nP Jf^{xK)J^[n)dn (LI) equals: 

• forpiTL, 

[ { (X - 1)1'' sin (.i-^) + {x- 1)7 sin (.i±^) } .FI (i-^, i-^; -p + 1; 1 - x'^) 

- sin (7ri±£^) \ I ^ ' \ ^ ^ ' x-^-P-^ 2FI (i±4^, ^+V^" ;p + 1; 1 - x'^) 

• forp&I. and i Z, 

jsin (7ri±^) [{X - + {-m^ - 1)7] + cos (7ri±^) (-1)'' 

PC l+p+ti+v l+p-p+v \ f 

+ xyS^.A^ - ^-') + (- 1)^ p| r-,^.+. j^| T3^j^ ^- ^-^-^ I sin (7ri±^) r„.„(l - x'^) 

- [sin (7r ^+^7+" ) logx-^(a; + 1) |x - 1| + cos (7r ^+^7+" ) n0{x - 1)] ai^i (i±£±i£±il, ^±£^±^1; p + 1; 1 - x'^) | 

w/zere aZZ f/ze equalities hold in the sense of distributions on R_|_. In both cases, the RHS is the sum of a 
distribution multiplied by a smooth function (the first term) and of a locally integrable function (the second 
term). We have denoted 9{x) the Heaviside theta function, and the functions 5'^,iy,p, T'/i.i/.p are defined for 
\z\ <lby 

a V 2 }k\ 2 ^fc fc 



2^ 

TT 



X a; 



00 n+p+p+iy \ j l+p-/j+l/ \ 

X + 1) + V(p + + 1) - V( '+^7+" +k)- ^( '+^7+" + fc)} , 
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with ip{y) := ■^T{y)/T{y). We use also the definition, for Re A > — 1; 




X < 1 



a; > 1 



X < 1 
X > 1 ' 



and extend it analytically in the sense of distributions to all values ofX€C (see Remark 4.2). 



2 The Weber-Schafheitlin integral as an integral kernel 

In the following section, we quote some results obtained in [BDG] for the Aharonov-Bohm hamiltonian and 
conclude from them, that the Weber-Schafheitlin integral describes the integral kernel of a physically relevant 
operator. We motivate thus the need for explicit formulae, vaUd in particular for the distributional Re p > 1 

case. 

We consider the Hilbert space L^(R^) and denote its inner product {■]■). Since it will be convenient to use 
polar coordinates r, ^ on K^, we introduce the unitary transformation 



L^{M.^) 3 f^Uf e 2.2(0, oo) ® L^i-TT, tt) 

given by Uf{r, (f)) = y/rf{r cos (j),r sin (j)), which allows us to identify ^^(M^) with ^^(0, oo) (Si L'^{—'k,tt). 
The Aharanov-Bohm hamiltonian in polar coordinates is 



understood as the self-adjoint operator associated to the differential expression above (defined on an appro- 
priate domain). We allow for the moment the parameter A to be any complex number. 

Since the self-adjoint operator L := —id^ has spectrum sp (L) = Z, we have the decomposition ^^(R^) = 

® Hk where Hk is the spectral subspace of L for the eigenvalue k. With the help of U we can identify Hk 

feez 

with L^{R). Since L commutes with H^^, we obtain the decomposition 



We now assume /z > — 1. We gather some results from [BDG] about the operator iJ^. We will need first 
to define the following synnmetric operator, corresponding up to a constant factor to the so-called Hankel 
transformation: 

Definition 2.1 is the operator on L'^{0, oo) given by 



2_ 1_ 

where H^j, acts as the differential operator —d'^+ ^i* , when restricted to Cc(M+). 




We have then: 
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Theorem 2.2 Let < a < b < oo and denote '^[a,b] the characteristic function of the corresponding interval. 

The integral kernel o/]l[a,fe] (-ff/*) 

'^[a,b]{H^,){x,y) = I ^/xyJfj,{xK)J^{yK)KdK, 
considered as a quadratic form on C^(R+), that is, explicitly: 







(/I (if^)/) = /_ l[„,,](K2)|(J-^/)(«)|^d« 

for any f G C^(M+). We may thus identify 
where Q is the self-adjoint position operator, and in consequence. 

Note that in particular, T^j, is a unitary involution. For any 7 e C, one gets 

y^^'p,-^ fj, ~ Q ^ 1 
and in the sense above, the integral kernel of is 

/•oo 

H^{x,y)= / K^'' ^/xyJ^{xK)J^{yK)KdK, 
Jo 

which can be expressed in terms of the Weber-Schafheitlin integral with p, = u. 

We quote also the following result concerning the wave operators for ff^, assuming now /x e M: 

Theorem 2.3 For p,,u > 1, the M0ller wave operators associated to H/j,, exist and 



t— >±oo 



In particular, the integral kernel of the operator Clj^^^H^ = H^fl'^^^ may be useful in calculations. We have 
by the above considerations: 

and its integral kernel is equal to 

Jo \ y Jo \y J 

where the last integral is of Weber-Schafheithn type with exponent 27 + 1 and argument | . 

3 The Weber-Schafheitlin integral with Re p < 1 

Proceeding as in [DF], we quote the following classic result [W] for the integral involving the modified Bessel 
function of the first and second kind, 1^ and K^: 



Lemma 3.1 For Re z > 0, \z\ > 1, Re (i^ + p + 1) > |Re one has 

JO 

Proof. The conditions for n, v, p, for the convergence of the integral, are established using asymptotic series 
for the Bessel functions of the corresponding type. 

For the derivation of the integral, we first rescale the variable n by the factor z~^, expand Iv{k,/z) as a power 
series and then use, assuming \z\ > 1, 

J k^{k)k^-' = J J e-«'^(«2-i)A'-^K/3-idudK = 2'5-2r(^)r(^) , 

where we have substituted for Kf^ the corresponding integral representation, and computed the obtained ex- 
pression, integrating first with respect to k. It remains to compare the obtained series with the hypergeometric 
2 -Fi function on the RHS . □ 

We recall the relation between the Bessel function of the first kind and : 



Using 



POO POO 

/ KPK^{zK)h{K)dK = Z-P-^ KPK^{K)h{K/z)dK, 

Jo Jo 

POO poo 

/ K''Ki,{zK)J4K)dK = Z-P-^ / KPK^{K)J^{Klz)dK, 

Jo Jo 
we get as a straightforward corollary of Lemma 3.1: 

Corollary 3.2 For Rez>0, Re (i/ + p + 1) > |Re iJ.\, 

POO 

/ KPK^{zK)MK)dK = T{ )r( )2P-iz-i-p-'- ^Fl l±p^t+Ji- -z-^) 

Jo 

(3.1) 

Note that the superfluous assumption |z| > 1 has been eliminated by analytic continuation with respect to z, 
using the analytic! ty of the Gauss hypergeometric function 2-F1 on C \ [ 1, 00 [. 

We denote the Hankel function of the first and second kind, respectively — H^, H~ , and recall that they are 
related to as follows, for any y gC: 

ITT 

H-{y) = -^e-^K^{iy). 

^ ITT 

It follows that the integrals 



POO 

/ KPH^{xK)J^{K)dK 

Jo 
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for X G M+ are both the Umiting case of n'' K ^{zk) Ju{n)'^n for purely imaginary z. We use the results 
obtained in Corollary 3.2 for Re z > 0, setting first z = ^i{x ± ie), which gives: 

POO yOO 

/ k,''H^{xk)Ji,{k)Ak = \\m I K^H^(x ±isK)J,^(K)dK 
Jo V y V y ^^^j^ ^ V 

= —62 lim / k''KJtHx ±ie)K)Ji,(K)dK 

ITT e\Oio 

= ±^e±^"'^'^r(i+£±ii+^^)r(i±£^)lim/(x±i£), 

ITT J, A 



where 



/(a; ± i£) := (a; ± xeY^-P-'' i±£fi±^^; + 1; (a; ± ie)-^) . 

Since we have also the relation = i + [W], it follows that 



/ K''j^ixK)j,iK)dK = ^r( ^+p+^+" )r( ^+^7+'- ) 



(3.2) 



Therefore, in order to derive the Weber-Schafheitlin integral, as well as the integrals involving H^, it is 
enough to examine the hmit limgN^o /(a^ ± ie)- Note that f{z) depends on the parameters /x, p and it will 
follow that it is convenient to treat the cases Re p < and Re p > separately. 

Proposition 3.3 For any /i, e C and Re p < satisfying Re {p + p + 1) > and x £ M+, the integral 

n'^ J ^{xn)Ji,{n)dn is equal to: 



for X < 1, and 



for X > 1. 



p/ l+p+H+i' \ 



Proof. Consider the factor 2^1 ( 1+^+^+^^ , i+p^^+t^ .^ _^ ^. _^ ie)-^)^ appearing in the definition of 
f{x+ie). Denotinga := ^+^+^+'' ,6 := andc := we checkthatRe (c-a-fo) = -Rep < 0, 

which ensures the limit £ \ exists. 

For a; > 1, the argument of the 2F1 function has real part smaller than 1, thus by analyticity it follows that 
the limits with +\e and — ie coincide. Obtaining the desired expression is then just a matter of rewriting the 
phase factors in terms of F functions, using 



1 ("e-^i^ - e--^^^"! = sin = 1 

2i V ) ^ 2 ) p ^ i+p+ix-M ^ r ( ^"P"'' 

For a; < 1, we can use the result above with fi and u interchanged, thanks to 

/•oo POC 

I Jij,{xK)J^{K)dK = X~''~^ / k'' J^(k) Ji/(K/x)dK. 

Jo Jo 



□ 
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4 The case Re ^ > 

Assuming Re p > 0, we examine the limit 

lim f(x ± i£) = lim (x ± iey^-P-" l+P^n+iL.y + i-(x± ie)-2) . 

Denoting a := i±£±li±il, 6 — and c := + 1, we check that Re (c - a - 6) = -Re p > 0, which 

implies a singular behaviour of the 2-F1 factor at a; = 1 in the limit £ ^ 0. We therefore use the following 
symmetry of the 2^^! function: 

-.Fl (a, 6; c; ^"2) = (1 - ^-2)c-a-b ^i^i (c - a, c - 6; c; ^"2) 
= (1 - ^-2)-" ( ^-^+^+" , IV + 1; ^-2) , 

in order to isolate the whole singular factor (1 — the second term on the RHS being well defined for 

all 2: e C. Consequently, 

j{x ± i£) = {x± i£)-i-''-^(l -{x± i£)-2)-^ -.Fl 1/ + 1; (a; ± ie)-^) 

= {x± xe)-^^P-''{x + 1 ± ie)-''(a; - 1 ± ie)"^ 2-F1' '"V^" ; + 1; (a; ± ie)"') . 

Because of the singular part (i.e. (a; — 1 ± ie)"''), the limit £ \ requires a distributional approach. Let us 
concentrate now on the first factors of the above expression. We will need the following: 

Definition 4.1 Wfe define {x— 1 ± 10)^ lim£\^o(2; — 1 ± ie)'*', as distributions on R+. 

Remark 4.2 The distributions [x — 1 ± 10)'^ are well defined and can be represented in a more explicit way 
as follows ( see [H] for properties of the analogously defined distributions {x ± 10)^ on Mj." 

• for Re A > —1, we have the equality between locally integrable functions 

(a; - 1 ± iO)^ = e=^'^''(a: - 1)^ + {x - 1)^, 

where: 

{x-lf_ :-- 



|a:;-l|^, x<l 
0, a; > 1 



0, x<l 
{x-l)^, x>l' 



for Re A > —k and A ^ Z, where A; e N, we have the equality 

(x - 1 ± iO)^ = e^'^^'ix - 1)^ + {x- 

where (x — 1)^ are now distributions defined by their action on an arbitrary test function (f> ^ 
(or equivalently, as the analytic continuation of the distributions defined in the preceding case): 

{{x - lf_,<l>) := /V - ((A + 1) ... (A + k)) dx, 

Jo 

{{x - 1)1, <t>) := (-1)'= {x - 1)^+^^') {x)/ ((A + 1) ... (A + k)) dx. 
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for \ = —k, where fc e N, we have 

{x-l± iO)-'= := {X - 1)Z' + {-l)\x - 1);'= ± {-l)Hn^^^j^., 

where 



((I -l):',.^> :=(-!)'-' 



i log(I-.)|;3^dx+(-l)'- L. 



(fc-1)! 

Lemma 4.3 Let A € C anc/ let g{z) be a function holomorphic in the neighborhood of the halfiine M+. Then 

lim \g{x + \e){x - 1 ± k)-^] = S((a;)(x - 1 ± iO)^, 

in the sense of distributions on 1R+. 

Proof. It is clear that 

lim [g{x + \e){x - 1 ± ie)^] = g{x){x - 1 ± iO)^, (4.1) 

for Re A > —1, since the RHS is then a regular distribution and the convergence of each of the factors as 
functions is uniform on every compact subset of M+. Assuming A 7^ 0, we differentiate (4.1) and obtain 

lim {g'{x + ie)(a; - 1 ± \e)^ + \g{x + i£)(a; - 1 ± \ef-~^\ = g'{x){x - 1 ± iO)^ + Xg{x){x - 1 ± iO)^-\ 

where we have used the fact that ^(a; - 1 ± 10)'*' = \{x - 1 ± 10)'*'"^ [H]. We then use (4.1) again to 
substract the first term of both sides. We have thus proved (4.1) for Re A > —2, A ^ —1. 

We consider the case A = — 1 separately, and differentiate instead the equaUty between the locally integrable 
functions: 

lim \g{x + \e) log(a; — 1 ± ie)] = g{x) \o%{x — 1 ± 10), 

where log(a; — 1 ± iO) := lim£\^o log(x — 1 ± ie) = log \x — 1 1 ± YnQ{x — 1), and its distributional derivative 
is (a; — 1 ± iO)~^ (this can be seen by setting first g(x) = \ in the above equality and differentiating). 

By induction, we prove (4.1) for the remaining values of A € C. □ 

Recalhng our expression for f{x + ie), we have 

f{x + ie) = {x + \s)-'^+P-''il +X + ieyix - 1 + ie)-Pq{x + ie), 

where 

q{x + is) := 2FI (i^£±^, +l;{x + ie)-^) , 

which converges uniformly on every compact subset of to 

q{x + iO) := ^-P-^+- -u + l;x-'), 

considered as a single-valued function. Note that the Gauss hypergeometric function 2F1 has a branch cut 
[ 1,00 [. Here, the limit is taken by approaching the halfiine from below (Im (cc + ie)~^ < 0), and we 
have denoted 2-^1 (a, 6; c; x) = limes^o 2Fi{(i. Ir. c; x — ie) on the branch cut. On the other hand, the limit 
limgx^o f{x — ie) corresponds to approaching the halfiine from above in the argument of the 2F1 factor, 
therefore 



q{x — 10) := lim q{x — ie) 



is not equal to g(a; + iO). 
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Remark 4.4 The function x i~> q{x ± iO) is not dijferentiable, and therefore the meaning of the product 
{x — 1 ± \Q)~''q{x ± iO) is unclear. To prove that it exists despite this apparent problem, we show that 
q{x ± iO) can be written as 

q{x±iQ) = hi{x) + {x-l±iQYK^ {x) , (4.2) 
where hi{x) is smooth and both belong to I/'°'^(M_|_). Then, the equality 

(a; - 1 ± iO)-'' q{x ± iO) := (a; - 1 ± iO)-^/ii(a;) + hfix) 

defines the desired product well, being the sum of a distribution multiplied by a smooth function and of a 
locally integrable function. 

Proof. For p ^ Z, the decomposition (4.2) is possible due to the following formula, holding for ^ e C ([BS], 
eq. (B.9)): 



2-^1 V 2 ' 2 



■v + l;z) = 1^ ^ ^ ^ 1^ ^ 2^1 ( '-Pt^+" , ■ p+i-i-z) 

^ sinTrp \r( ^+^+^"'"'' )r( ^+^~^+'' ) ^ ^ ^ > p I > ; 

- (1 - 0)^ 2FI ( 1+^+^+-- , 1+^-/^+'- .p+i-i-z) |.(4.3) 



We take z = {x± ie)"^ and pass to the limit £ \ with both the expressions. Note that the only difference 
between the limit with +ie and —\£ appears in the factor lime\^o(l — (a;±ie)~^)'' = a;~^^(l+x)''(x— l±iO)'', 
since the 2^1 factors on the RHS are analytic on the required domain. It is clear that (4.2) holds, with 

^2 (^) = ^2 (2^) ill particular. 

For p e Z, we have instead ([BS], eq. (B.IO)): 

1 ^ {^Ifip -k-l)\ (^^^), (^"^), 

_^(i±p+i£±ii + _ ^( 1+^7+- + - iog(i - z)] (1 - 

for i^£±if±;i ^ z, where ^(j/) = T'{y)/T{y). 

As in the preceding case, we take z = (x ± ie)~^ and pass to the limit. We obtain a similar decomposition. 
Note, however, that the limits log(l — .t =F 10) differ, and thus h^{x) 7^ h2 (x) 

The cases where at least one of the parameters ^~p^t^+'' is an integer are treated the same way, only the 
functions hi{x) hf{x) being then different ([BS], eq. (B.ll), (B.12)). □ 

Proposition 4.5 In the sense of distributions on M.+, 

lim /(x ± i£) = x-'^+P-''{x + 1)-P [(a; - 1 ± iO)-Phi{x) + hf{x)] , 

where the functions hi{x) and /i^(x) are defined in Remark 4.4. 

Proof. We have to prove that as £ \ 0, 

(x - 1 ± i£)-''(x ± i£)-^+''-''(x + 1 ± \e)-Phi{x ± i£) ^ (x - 1 ± iO)-''x-^+''-''(x + l)-''/ii(x) 
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and 

{x ± \e)-^+P-''{x + 1 ± leyix - 1 ± \e)-P{x - 1 ± iefhiix ± k) — )• a;-^+^-''(l + x)-Ph^{x). 

The first limit is a consequence of Lemma 4.3. The second one is clearly true, since the convergence of the 
corresponding functions is uniform on each compact subset of R+ . □ 



Remark 4.6 The case when p e Z and at least one of the numbers is an integer are treated 

similarly. One can deduce from the expansions given in [BS] (eq. (B.ll), (B.12)), the explicit expressions for 
the functions h\{x) and hf{x), following step by step the proof of Remark 4.4 in the degenerate case. 

Corollary 4.7 Using Equation 3.2 and Proposition 4.5, we get the results gathered in Proposition 1.1. The 
formulae in Proposition 1.1 are written in terms of the distributions {x — 1)^^ rather than (a; — 1 ± iO)~^, 
using the relations listed in Remark 4.2. 

We end up commenting on some special cases, involving much simpher expressions. 

Remark 4.8 An explicit formula in the special case p = I has been derived in [KR]. It can be recovered 
from our general expression, by substituting the (well-known) equality 

(a; - 1 ± iO)-i = Pv {^-^—^ T i7r(5(a; - 1), 

where Pv denotes the Cauchy principal value. Furthermore, S'^,i/,i(a;) = 1 by definition, and it remains to 
use (4.3) back again to get 2F1 functions with argument x"^^ instead of 1 — a;^^. 

Remark 4.9 Much simplier formulae can be derived in the special case 1 — p + p.±v = Q, since the Gauss 
hypergeometric function with a parameter set to zero is trivial, i.e. 2-F1 (0, . . . ; •) = 1. 

Note that the above remarks hold for the integral involving H"^ instead of as well. 
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